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Feedback in the Antennae Galaxies (NGC 4038/9): I. 
High- Resolution Infrared Spectroscopy of Winds from Super Star 

Clusters^'^ 
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ABSTRACT 

We present high-resolution {R ~ 24, 600) near-IR spectroscopy of the 
youngest super star clusters (SSCs) in the prototypical starburst merger, the An- 
tennae Galaxies. These SSCs are young (3 — 7 Myr old) and massive (10^ — 10"^ 
Mq for a Kroupa IMF) and their spectra are characterized by broad, extended 
Brackett7 emission, so we refer to them as emission-line clusters (ELCs) to distin- 
guish them from older SSCs. The Bvy lines of most ELCs have supersonic widths 
(60 — 110 km s~^ FWHM) and non-Gaussian wings whose velocities exceed the 
clusters' escape velocities. This high-velocity unbound gas is flowing out in winds 
that are powered by the clusters' massive O and W-R stars over the course of 
at least several crossing times. The large sizes of some ELCs relative to those of 
older SSCs may be due to expansion caused by these outflows; many of the ELCs 
may not survive as bound stellar systems, but rather dissipate rapidly into the 
field population. The observed tendency of older ELCs to be more compact than 
young ones is consistent with the preferential survival of the most concentrated 
clusters at a given age. 
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Suhject headings: galaxies: individual (NGC4038/9, Antennae Galaxies) — galax- 
ies: ISM — galaxies: starburst — galaxies: star clusters — infrared: galaxies — 
HI I regions 



Introduction 



Star formation in starbursts creates massive young super star clusters (SSCs) that are 
not often found in more quiescent environments like the Milky Way's disk. Their inferred 
masses (10^ — 10^ M©) and Lyman continuum photon rates (10^^ — 10^^ s~^) dwarf those of 
their lower-mass analogs in the Local Group. The most massive local SSC analog is R136 
i n 30 Doradus, with a mass of at least 10^'^ M© (IMassey fc Hunterl Il998al ) and 10 
flWalbornlll99lh . 



51.4 g-1 



NGC 



1995 



1995 ) , to merging sy stems such as the Antennae Galaxies (NGC 4038/39) (e.g. lWhitmore fc Schweizer 



. O'Connell et al. 




Whitmore & Schweizer 



Larsen 



GilbertI 120021) . and some have been found in normal spiral galaxies (like NGC 6946 
Larsen et al. boOl), wh ose young clusters typically populate a lower mass range, ~ 10^ — 10^ 



Mq (e.g. 

but near-infrared (near-IR) imaging has revealed young massive clusters t hat excite lu- 



2OO2I ). Most observations of SSCs have focused on optical wavelengths, 



minous H II regions (H IIRs) that are invisible at optical wavelengths (e.g. [Gilbert et al. 
2OOOI : [Turner et al.ll2003l ). Narrow-band imaging may also be important to track the faintest 
youngest members of the population, which may not be detected with broad-band IR imag- 
ing although they are bright in narrow filters tuned to the recombination lines of their HI- 
IRs (lAlonso-Herrero et al.ll2002l ). Radio and mid-IR observations have uncovered still more 



heav ily embedded regions of mass i ve star forrnation that are heavily obscured even at 2 /xm 



(e.g. 


Kobulnickv & Johnson 


1999: 


Beck et al. 


2000. 


2001; 


Vacca et al. 


2002: 


Cabanac et al. 


2005 


); these objects are very young, dense, compact HIIRs that are the massive analogs of 



ult ra-compact H ll regions (UCHIIs). Hence they are dubbed ultra-dense HIIRs, or UDHIIs, 



by iKobulnicky fc JohnsonI (119991 ) , and they are the supposed precursors of optically visible 
SSCs. 

The high masses and stellar densities of SSCs resemble those of globular clusters (GCs), 
although metal abundances ought to be higher in SSCs than GCs and SSC ages are mea- 
sured in Myr rather than Gyr. This suggests an evol utionary sequence in massive star- 
forming regions: UDHIIs become optically visible SSCs (IKobulnicky fc JohnsonI Il999l ) that 
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may finally become (metal-rich) GCs. The idea of GC formation in mergers is su pported by 



obseryations of young GCs ( ages of a few Gyr) in several merger renmants (e.g. ISchweizer 



1987 



Ashman &: Zepj Il992l : iFritze-v. Alvensleben fc Gerhard! 1 1994 ISchweizer et al 



1996 



Schweizer &: Seitzer 19981 ). as well as the multi-modal color distributions and varied spatial 



distributions of GCs. However, estabhshing the relationship between SSCs and GCs requires 
examination of both the photometric and kinematic properties of the SSC population. By 
measuring an individual SSCs age, mass, density profile and initial mass function (IMF), and 
modeling its stellar and dynamical evolution, one can predict whether it may survive internal 
and galactic dynamical processes to 



Boker et al 



2002 



1999 



oin a population of GCs (e.g. IHo &: Fi 



Smith fc Gallagherl l200ll : iLarsen et al.l 1200 ll : iMengel et al 



ipperiko 



2002 



1996 



Gilbert 



McCrady et al.ll2003l ). The IMF of SSCs is especially important for this question be- 



cause the light of SSCs is dominated by massive stars, while old GCs are comprised entirely of 
sub-solar mass stars (and stellar remnants). Observations of SSC dynamical masses to date 
sugg est that some SSCs do have normal IMFs and could become GCs, while others cannot 



McCrady et al. 



2003 



Sternberd 19981: ISmith k Gallaghei]l200ll : iMengel et al.ll2002l : ICilbert k Grahamll2001 



200 



^, although these results could be subject to systematic errors in 



extinction corrections or neglect of effects like spatial mass segregation, and they always rely 
on the assumption of cluster virialisation. 

Another important issue is the evolution of the SSC population as whole: GC popula- 
tions have kinematic distributions resembling those of bulges (spheroidal with high velocity 
dispersions), while SSCs form where gas is available, usually in dis k-like enviro nments with 
much lower dispersions; GCs have log-norm al luminosity functions (jHarris 1991), while SSCs 
display power- law luminosity functions (e.g. lWhitmore et al.lll999l ). as do the giant molecular 
clouds (GMCs) and H ll regions from which they pres umably form (e.g. IWilliams k McKee 
19971 : iKennicutt k Chulll988l : iMcKee k Williamslll997l ) . Simulations have shown that power- 
law cluster LFs can evolve into log-normal ones as the faint end of the LF is eroded by tidal 
disruption in the galacti c potential, dynami cal friction, stellar mass loss, and evaporation of 



stars from clusters (e.g. iFall k Zhang l200ll ). 



Whether young SSCs survive to become GCs or disperse completely into the field star 
population of a galaxy, they have great potential to affect the energetics of its interstel- 
lar medium (ISM) because they harbor thousands of massive stars producing ionizing and 
FUV radiation. The radiation from OB stars excites HIIRs and photodissociation regions 
(PDRs), and heats the ISM. Their winds and supernova (SN) ejecta stir and inject energy 
into the surrounding ISM. In the most extreme st arbursts, the combined effects of a star- 
burst drive large-scale galactic winds (e.g. in M82. IShopbell k Bland-Hawthorrull998l ) that 
can eject matter into the intergalactic medium (especially in dwarf galaxies, e.g. NGC 1569, 
Heckman et al.lll995l ). SSCs are concentrated power sources for feedback on both star-cluster 
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and galactic scales. 



The Antennae (NGC 4038/9) are a nearby pair of disk galaxies in an early stage of 
merging that are well-known for th eir numerous SSCs, d istributed along their spiral arms 
and around their interaction region Jwhitmore et al.lll999h. Their mole cular gas distribution 
peaks at both nuclei and in the overlap region (IStanford et al.lll990l ), but the gas is no t 
yet undergoing a global starburst typical of more advanced mergers (INikola et al.lll998l ). 
The current global star-formation rate of the system is estimated at about 20 Mq yr~^ by 
Zhang et al.l (120011 ) from extinction-corrected Ha. ISO observations of the Antennae show 
that the mid-IR (8 — 15 /im) flux from warm dust follows the optical distribution of bright 
blue SSCs, but peaks in t he overlap region in a single point source that emits 15% of the ISO 



flux (Mirabel et al 



19981). T his source is an optically faint and unremarkable star cluster 
that near-IR observations reveal as one of the voungest f~ 4 



(Gilbert et al. 


2000; 


Mengel et al. 


2001) 



Myr), most massive (~ lO'^ M©) SSCs in the systi 
This cluster powers a dense (average Ue = 10^), large (half-light radius 32 pc) H ll region 
with a Lyman continuum rate of order Q[H^] = 10^^ s~^; it excites even more extended, 
clumpy PDRs with nearly fluorescent H2 emission; and it is embedded behind Ay ~ 10 mag 



(IGilbert et al.ll2000l ). 



In § [2] we present near-IR imaging and spectroscopic observations of a sample of young 
Antennae SSCs in order to identify the emission-line clusters (ELCs) among them. We derive 
ages and stellar masses for the ELCs from their Br7 spectra and magnitudes in § [31 In § H] 
we compare the widths of the broad, extended Br7 emission lines of ELCs with the inferred 
escape velocities of the exciting SSCs, and argue that ELCs drive cluster outflows. In § [5] 
we estimate electron densities for ELCs, and we compare ELCs with other types of HIIR. 
We discuss the evolution of ELCs and survival of massive clusters undergoing expansion due 
to mass loss in § El estimate the ELC contribution to the total star-fo rmation rate in the 
Ante nnae in § [71 and conclude in § [HI In a companion paper (Paper II, [Gilbert fc Graham 
20071 ). we present a kinematic model for the cluster outflows, fit it to the data presented 
here, and discuss the mass-loss rates and energy budget of ELC outflows relative to similar 
systems such as superbubbles and galactic superwinds. 



Observations &; Data Reduction 



On 2002 January 21 and February 22, we obtained i^'-band spectra of a sample of 17 
young SSCs arid the nuclei of NGC 4038 and NGC 4039 in the Antennae using NIRSPEC 
([McLean et al.l [20001). a facilit y near-infrared (0.95 — 5.6/im) spectrometer for the Keck-II 
telescope ([McLean et al.[[l998[ ). Using the cross-dispersed echelle mode and the N7 order- 
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sorting filter, we detected Brackett7 with little or no continuum in 16 of the targets with 
A/ A A ~ 24,600. During the course of these observations we recorded N7 images of the 
fields around the 0."432 (3-pixel) slit using NIRSPEC's slit- viewing camera (SCAM); these 
images were used to create a mosaicked finding chart of the Antennae field from which offsets 
between the infrared clusters were measured. Cluster photometry was performed on iiT-band 
images taken with NIRSPEC's slit-viewing camera (SCAM) in good seeing (0."6 or better) 
and photometric conditions on 2000 January 17 a nd 2001 March 10. We later obtained 
higher quality VLT ISAAC jMoorwood et al I llQQsh narrow-band images at 2.17 and 2.25 
/im on 11 January 2004 in K-band seeing of 0'.'4 — O'.'S; the former filter includes Br7 emission 
and is shown in the finding chart of Figure [H 



2.1. Imaging 



The bright northern and southern sources in Figure [T] are the two nuclei in the system, 
NGC 4038 and NGC 4039. Overplotted are the positions of the brightest young (< 30 Myr 
old) clusters detected in optical HST images reported by lWhitmore et al.l (119991 ). The HST 
and NIRSPEC data sets were registered using the foreground star at the origin of the figure 
(Star 4 of Whitmore et al. 1999)lj. Comparing the relative HST positions of the nine cleanly 
overlapping sources with their relative IR positions yielded a S CAM plate scale of 0'.'178 
per pixel. This is consistent with the value of O'.'IS measured by iFiger et al.l (120001 ). About 
half of our sample of bright young SSCs does not appear in the optical brightest cluster list 
due to dust extinction, although most of these sources are detected at least in F814W (see 
Figure [2]). SSC Bl highlights the importance of long- wavelength observations in probing 
the most recent star formation: it is deeply embedded an d quite faint in the optical, with a 
K-band extinction of at least 1 mag (IGilbert et al.ll2000l ). and it is the most massive young 
SSC in the system after correcting for this extinction (Table H]). 



The i^'-band SCAM images were reduced and calibrated relative to the iPersson et al. 



( Il998l ) standards SJ 9150 and SJ 1952. Aperture magnitudes for the clusters were calculated 
using DAOPhot (IStetsonlll987l ) routines in IDL. The fields containing the clusters are complex 
and it is evident that the accuracy of our photometry is limited by systematics associated with 
eliminating contributions from unrelated sources and with estimating the local background. 
In an attempt to reduce these systematic errors we also computed magnitudes based on 



^ Whitmore et all (|l999l ) report a position for Star 4 of RA = 12''1'"56.04'' and DEC = -18°52'43'.'66 
(J2000), althou g h an offset of 1'.'2 t o the southwest of the H ST positions from radio observations 



Neff k UlvestadI (|2000l) was reported bv lWhitmore fc Zhanel (|2002f) 
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fitting a 2-d Gaussian to the resolved cluster light profile. The mean difference between the 
2" diameter aperture and fitting magnitudes is 0.12 mag with an rms of 0.2 mag. Inspection 
of the field confirms that the 2" apertures represent a compromise between isolating the 
clusters from their environment while enc ompassing most of the light (point-source aperture 
corrections were applied using Star 4 of IWhitmore et al.l 119991 ) . We therefore adopt these 
magnitudes in what follows and assume 0.2 mag as the rms error. Table[T]lists the magnitudes 
of the selected targets for spectroscopy. Our photor netry agrees with o ther measurements 
from the ht erature: for SSC Bl (known as SSC A bv Icilbert et~aD boooh . w efind K = 14.7 
± 0.2 mag, ICilbert et all kood ) report N7 = 14.6 mag, and iMengel et all J200lh r eport K 
= 14.8 ± 0.2 mag; for SSC F we find K = 15.9± 0.2 mag and iMengel et all J2002h report 
Ks = 15.9 n iag. For SSC S we d id not have adequate data for photometry, so we use values 
reported by IMengel et al.l (120021 ) . 



The ISAAC images were reduced, corrected for distortion, and combined. Sizes of ELCs 
and Star 4 were measured via 2D Gaussian and Lorentzian fits, which gave comparable 
results, and the latter were adopted because they better reproduced the observed profiles. 
We derived FWHMs from the geometric mean radii of the fits because it is a better indicator 
than the arithmetic mean of the area in an image that encloses half of the profile's light. The 
FWHM of Star 4 was subtracted in quadrature to compute the deconvolved FWHMs that 
are reported in Table [1] and plotted in Figure [61 Throughout this paper we adopt a redshift 
distance for the Antennae of 19.2 Mpc (1" = 93 pc, distance modulus (m — M)o = 31.41 
mag for Hq = 75 k m s~^ Mpc~^), although we note that a smaller distance of 13.8 Mpc has 
been measured by ISaviane et al.l (|2004l ) based on the tip of the red giant branch (TRGB) 
in the southern tail of the Antennae. If the TRGB distance is correct, all linear dimensions 
discussed in this paper decrease by a factor of 1.4. 



2.2. Brackett7 Spectroscopy 



A finding chart similar to that in Figured] from prior i^'-band observations permitted the 
selection of cluster targets for spectroscopy. Targets were chosen on the basis of brightness; 
whenever feasible we oriented the 24" long slit to encompass multiple clusters. 

Table [2] lists exposure times for the targets that are identified in the image of the 
Antennae in Figure [H SSCs are named ro ughly accordi n g to t he lettering scheme adopted 
for Ha-emitting regions in the system by iRubin et al.l (|l970l ). Along with the SSCs we 
observed the two nuclei and detected broad emission in the northern nucleus, NGC 4038, 
but not in the southern one, NGC 4039, although we did detect broad Br7 in a faint cluster 
near the latter nucleus. This off-nuclear source is labeled "Al" in Figure [Hand in Table [31 
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Raw exposures were first corrected for electronic noise which appears as horizontal 
striping in a few of the 32 analog channels. The images were then dark-subtracted (because 
on-minus-off sky subtraction was not possible for a few exposures) , fiat-fielded, and corrected 
for b ad pixels and cosmic r ays. A wavelength scale was determined from atmospheric OH 



lines (IRousselot et al.ll2000l ) and used to rectify the Brackett 7 echelle order (order 35), in 
which spatial and spectral directions are curved, onto an orthogonal grid of wavelength and 
slit position. In order 35 the spectral resolution measured from OH lines was X/6X = 24, 600, 
or about 12 km s'^ FWHM. 

A primary sky subtraction was performed by differencing pairs of frames where possible, 
and then a secondary sky subtraction was done by fitting and removing the background in 
each column. Emission-line spectra were optimally extracted using a Gaussian weighting 
function matched to the spatial extent of the Br7 emission, and then an aperture correction 
was applied to recover the full fiux. 

A B4IV/V star was observed as an atmospheric standard. Although most of our Br7 
spectra were unaffected by atmospheric absorption, at the radial velocities of some objects 
a small feature absorbs some flux in the wing of the line, so we applied the atmospheric 
correction in all spectra. A flux scale was derived by requiring the 2.2 fim continuum flux of 
the standard star to equal that corresponding to its K magnitude. Reduced line spectra for 
the ELCs from which we detected broad Br7 (excepting only Dl) are shown in Figure [3l and 
position- velocity diagrams for the four brightest ELCs are shown in Figure HI Signal-to-noise 
ratios per pixel in the fully reduced spectra ranged from around 3.5 in the weakest lines to 
75 in the strongest, with continuum levels low in all but the brightest sources, which have 
continuum signal-to-noise ratios per pixel of only a few. Sources for which we did not detect 



Br7, clusters 10, 16, and 1 from IWhitmore et al.l (119991 ). are not considered further in this 



paper. We refer to the SSCs with Br7 emission as emission-line clusters (ELCs). 

The Br7 proflles in Figure [3] are well resolved and exhibit high-velocity non-Gaussian 
wings. The isothermal sound speed of ionized gas at temperature T = 10^ K is c = 
ksT / ^ 12 km s~^ where /i is the mean molecular weight and mp is the proton mass, 
so linewidths FWHM = 2.355c > 28 km s~^ are supersonic. Nearly all of the line widths, 
determined from Gaussian flts to the line cores, are supersonic. Cluster Dl (FWHM = 14 
km s~^) is the only Br7 source without a supersonic component and cluster F is marginal 
with FWHM = 30 km s~^ We refer to these young SSCs with supersonic Br7 line widths 
as broad emission-line clusters (bELCs). 

The broad Br7 lines of ELCs may arise from individual hot stars or dilute p hotoionized 



gas. Individual O stars have EW[Br7] of only a few A (e.g. iHanson et al.lll996 ). while the 



Antennae clusters have EW[Br7] from a few up to a few hundred A. The observed Br7 lines 
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are also too narrow to be associated with individual 0-star winds, and as already noted, 
they are too broad (and non-Gaussian) to simply reflect the virial motion of the constituent 
stars. Thus, photospheric and wind emission from O stars alone cannot explain the intensity 
or width of Bvy emission from ELCs. Even though the ELC lines are too n arrow to be 



stellar wind features, WR stars could produce the observed equivalent widths. iFiger et al. 



( 119971 ) measure the combined EW of blended He I, He ll and Br7 lines from low-resolution 
IR spectra of WR stars, finding values of about 20-130 A. However, the line widths are 
far too narrow to be dorn inated by WR emission (which has line widths ~ 10^ km s^^. 



e.g. lAbbott fc Contil 119871 ). Thus the ELC line fluxes most likely originate not from gas in 
photospheres and winds of hot stars, but in the expanding H ll regions that they excite. The 
broad lines must be formed from rapidly expanding ionized gas swept up by the cumulative 
effect of stellar winds and photoevaporated flows driven by the thousands of hot stars in the 
clusters. 

Table [3] summarizes the measured Br7 fluxes, equivalent widths, barycentric radial 
velocities, and line widths for the ELCs. Our Br7 radial v elocit ies are consistent with those 
reported from HST Ha spectroscopy by IWhitmore et al.l (120051 ) for sources in the vicinities 
of SSCs B, D, and F (a direct cluster-to-cluster comparison is not possible because the sht 
positions and sizes are different). The listed Br7 FWHMs are the observed values, not 
corrected for either the instrumental line-spread function (12 km s~^) or thermal broadening 
(21.4 km s~^). We also detected He I 3F°— 3G emission at 2.1647 fim near the blue wing of 
Br7 in the brightest sources (Bl, B, C, and D) at a level of 4.4 ±0.5% of the Br7 fluxes. For 
SSC D2, whose Br7 emission includes a very narrow component (FWHM = 23 km s~^) as 
well as a weak broad component (FWHM = 55 km s~^), we list the total line flux but only 
the line width of the broad part. 

We consider the kinematics of the ELC population in the Antennae by computing the 
rms radial velocity dispersion of the 14 targets in the overlap region. We find a value of 
74 km s~^ rms, which is comparable with a disk galaxy's rms integrated velocity dispersion 
( lBottemalll993l ). Some of this velocity dispersion is due to a ge neral gradient in the north- 
south direction along the overlap region between the two nuclei ( lAmram et al.lll992l : [Gilbert 
2OO2I ). but the gr adient is not monot onic: both nuclei have higher radial velocities of at 
least 1600 km s~^ ( Amram et al. 1992[). and some ELCs lie at minirna in Ha radial velo city, 
i.e. in Ha superbubbles (lAmram et al.lll992l : IWhitmore et al.lll999l : IZhang et al.ll200ll ). In 
the northern overlap region, the ELC rms radial velocity dispersion is 21 km s~^ (for eight 
ELCs) about a mean of 1573 km s^^. In the southern overlap region, the velocity dispersion 
of six y ounger ELCs is 48 k m s^^, equivalent to that of the giant molecular complexes, (49 



km s-\ IWilson et all 120001 ) about a mean of 1445 km s ^, so the kinematics of this group 



reflect those of the nearby molecular material. Velocity dispersions for smaller groupings on 
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scales of 10 — 20", e.g. the two B, three D, five E, or three F clusters, are 21, 56, 22, and 18 
km s~^, respectively, which apart from the larger D cl uster value are co r npara ble with the 
raw dispersions for sources over similar spatial scales in lWhitmore et al.l (120051 ). 



High-resolution near-IR spectroscopy has recently revealed broad s upersonic Br7 emis- 
sion from massive embedded star-forming regions in two other systems: iTurner et al.l (120031 ) 
observed a deeply embedded radio source in NGC 5253 and iHenry et al.l (120071 ) observed 



two similar sources in Henize 2 — 10, for which they derived Qp"*"] values above 10^^ s~^ 
and FWHMs of order 60 — 80 km s~^. While these dwarf galaxies are different in scale from 
the merging Antennae disks, their starbursts have much in common: optically visible SSCs 
that have already disrupted their ISM on large scales (seen e.g. in X-ray and Ha superbub- 
bles), bELCs that appear to be blowing out of their natal cocoons, and even younger heavily 
enshrouded star-forming regions seen as radio-bright UDHIIs. 



3. Ages & Masses of ELCs 

Br7 fiuxes trace the ionized gas within a cluster and hence the young, massive stellar 
content. We assume that this fiux arises from the 0/WR stars of a population that formed 
in a short-lived burst and has an initial mass function (IMF) that is usually characterized 
as a power la w, MdN/dM oc M~" . Employing a stellar population synthesis model such as 
Starburst99 (ILeitherer et al.lll999l ) permits an inference of ELC age by comparing observed 
and predicted Br7 equivalent widths (EWs) for a cluster as a function of age, and an inference 
of mass follows from scaling the predicted magnitude of the fiducial 10^ M© cluster at the 
determined age to match the observed magnitude. We assume solar metallicity for the 
cluster^ and adopt the Galactic field star average IMF of Kroupa (a broken power l aw with 



a 



2.3 for stellar masses M = 0.5 - 100 M© and a = 1.3 for M = 0.1 - 0.5 M©, iKroupa 



200ll ). The slope of this IMF at the high- mass end is nearly equival ent to that of t he solar 
neighborhood IMF of Salpeter (a single power law with a = 2.35. ISalpeterl Il955l ). and is 
consistent with that measured by iMassey &: Hunterl (Il998al ) for the massive, dense cluster 
R136 in the Large Magellanic Cloud. 

Those assumptions yield ELC ages of 3.4 — 7 Myr and masses of 7 x 10^ — 12 x 10^ 
Mq (incorporating an extinction correction where possible), as listed in Table HI The ELC 
ages have a gradient across the overlap region: the five ELCs that are younger than 6 Myr 
reside in the southern overlap region where the bulk of the molecular gas is concentrated 



We note that while iMengel et all |2002t ) found a handful of Antennae SSCs to be consistent with havine; 
solar metallicity, they argue that one (ELC S, their [W99]-2) may have a supersolar metallicity. 
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(e.g. IWilson et al.ll2000l ). and the older sources reside in the northern overlap region. We 
estimate fit errors for the ages of less than 0.1 Myr and relative errors in masses of 20% for 
most ELCs (see Table H]). However, more significant than the fit errors may be systematic 
errors, possible sources of which include the unknown i^-band extinction for some ELCs, 
our assumptions that the ELCs are ionization-bounded, and uncertainties in the models. 
We discuss these possibilities in turn below, and then compare our ages for two ELCs with 
literature values. 

From an ELC Br 7 line flux we derive QfH"*"] under case B assu mptions for an ionized H ll 



region at 10^ K (e.g. lOsterbrockl Il974j : iHummer &: Storeyl 119871 ). The resulting values are 



listed in Table H] where they are uncorrected for extinction, whereas in the text we always 
adopt the listed extinction corrections where available. ELC Bl has Q[H+] = 1.5 x 10^'^ 
s~^, so it emits nearly half of the total Lyman continuum rate from all 17 ELCs observed: 

3served for SSCs 
UDHIIs 



3.3 X 10^^ s ^. ELCs Bl and B have Q[H+] at the high end of the range o 



^52-53 g-1 



Johnson et al 



2000). 



Beck et all 120001 1. The faintest 



in other systems (e.g. He 2 — 10 region A at 10 
are also observed to have Q[H+] values of order io^^~^^ s~^ ( 
Antennae ELCs in our sample have Q[H"'"] = 10^^ s~^, which is comparable with the fluxes 
of the brightest Galactic giant H ll regions (e.g. W43 and NGC 3606). 

ELCs can suffer significant extinction even at K b and, and where estim ates are available 
we include them in Table HI In the case of ELC Bl, [Gilbert et al.l (120001 ) measured an ex- 
tinction at K band of = 1.1 ± 0.1 mag from Br7 and 20 Pfund recombination liri es. This 
corresponds to A^J = 10 mag for a standard extinctio n law (IRieke fc Lebofskvl 119851 ). Other 
authors have inferred lower extinctions t o this source: Mengel et al.l ( 20011 ) derived Ay = 4.3 



mag by comparing Br7 and Ha fluxes; IWhitmore fc Zhang (120021 ) report Ay = 7.62 mag 
from optical imaging. We use radio continuum flux densities for ELC counterparts (from 
Neff fc UlvestadI |2000| ) together with Br7 fluxes to estimate i^-band extinctions for a few 
ELCs (Bl, B, D2, and Fl) under the ass umption that the ra dio emission is unabsorbed. For 
ELC Bl, our inference of Ay = 10 mag ([Gilbert et al.ll2000l ) from K-band data agrees with 
the extinction derived using the radio-IR flux ratio, suggesting that the emitting nebular 
gas and obscuring dust in the ELC are probably mixed; optical measurements only probe 
part- way into th e emitting reg i on so they underestimate the overall extinction and hence the 
intrinsic (5[H+]. [Turner et al.l (120031 ) find the same bias for the compact, heavily obscured 
source in NGC 5253: near-IR and radio measurements give a consistent value for the extinc- 
tion of Ay = 18 mag, a value that is six times greater than the optically-determined one. 
We emphasize that infrared-radio determinations of Q[H+] are more reliable than optical- 
infrared ones. In Tables [3] and H] we present observed fluxes and derived masses (uncorrected 
for extinction). Extinction does not affect the ages since they are derived from extinction- 
independent equivalent widths. 
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There is evidence that ioniz ing photons escape from young clusters i n starburst HIIRs 
(e.g. in He 2-10 and NGC 4214.|johnson et aLlboOol: lleitherer et al.l |l996l) . Photoionization 



mode ls of GHIIRs require 10 — 73% of the Lyman continuum flux to escape (ICastellanos et al. 



2OO2I ). The escape frac tion of Lyraan co ntinuum photons from Galactic OB associations is 
estimated at 6 — 15% (IDove et al.ll2000l ). Observations of diff use Ho; emission throu ghout 
the Antennae (not only in association with stellar clusters, e.g. IWhitmore et al.lll999l ) may 
also sugges t leakage of Lyman co ntinuum flux from its HIIRs, i.e. that ELCs are density- 
bounded. iBeckman et al.l (120001 ) flnd that a transition between ionization-bounded and 
density-bounded HIIRs in spiral galaxies occurs at a constant Ha luminosity of L-^^ = 10^^-^ 
erg s~^. For Case B recombination this corresponds to a Br7 luminosity L-q^^ = 10^^-^ erg 
s~^, which is below the observed values for all of the ELCs. Presumably such a threshold in 
a merging pair of disk galaxies is higher than in a spiral because it depends on the available 
gas supply and the strength of the starburst, both of which are stronger in the overlap region 
of a merger. If ELCs are density-bounded, then some of their Lyman continuum photons 
escape to ionize and dissociate the surrounding medium, and our Br7 measurements would 
underestimate Q[H^]. We would thus overestimate ELC ages and expect them to be fainter, 
which would also lead to overestimated masses. 



For two sources in our sample, iMengel et al.l (120021 ) measured fC-band magnitudes and 
dynamical masses based on high-resolution spectroscopy, and used additional diagnostics 
with EW[Br7], i.e. EW[CO] and EW[CaT], to derive ages for the clusters. We compare our 
ages with theirs to give an idea of the significance of IMF choice, metallicity, and the level of 
systematic errors that may be inherent to the observations and Starburst99 models, which 
we both employ. For ELC S (their cluster [W99]2), we infer an age of 6.6 ± 0.1 Myr using 
EW[B r7], assuming a full (0.1 — 100 Mq) Kroupa IMF at solar metallicity. iMengel et al. 
J2002h obtain the same age (6.6 ± 0.3 Myr) for a Salpeter IMF (1 - 100 Mq) at twice solar 



metallicity. Our measurements of EW[Br7] differ by a factor of two (ours: 4.8 ±1.6 A, theirs: 
10.0 ±2.5 A), most likely because we use different apertures: theirs is a 2'.'2 box that includes 
strong diffuse line emission surrounding ELC S (which is not at a peak in Ha or Br7), while 
our aperture encompasses less diffuse line emission because it is the best-fit Gaussian spatial 
profile of the cluster Br7 line, or 0'.'9 FWHM for ELC S. In this case the differences in ages 
due to these other differences (between IMFs, abundances, and even EW[Br7]) cancels out 
to be insignificant. However, the typical size of the age offsets due to these factors is a few 
tenths of a Myr in the sense that the Salpeter IMF gives an older age than the Kroupa IMF 
for a given EW[Br7], and a larger EW[Br7] for the same IMF gives a youn ger age, while 



the rn etallicity difference in age is quite small for this range of EW[Br7]. But IMengel et al. 



( I2OO2I ) infer that neither of these IMFs may be appropriate for the cluster based on the 
dynamical mass measurement (2.0 x 10^ Mq), which is lower than the photometric mass 
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predicted by the models. This imphes t hat the cluster has less than a full complement of 
low-mass stars or a shallower IMF slope (IMengel et al.ll2002l ). 



In ELC F (their [W99]15), Mengel et all (120021 ) measure a high EW[CO 2.29/im] of 
17 ± 0.2 A. This requires an age > 7 Myr, when the brightest supergiants dominate a 
cluster's infrared emission. They estimate an age for ELC F of 8.7 ± 0.3 Myr for a Salpeter 
IMF, while we find 6.5 ± 0.1 Myr for a Kroupa IMF or (6.6 ±0.1 Myr for the same Salpeter 
IMF), where we both assume solar metallicity. Again the dynamical mass implies that 



neithe r of these IMFs is correct, but in the opposite sense to that of ELC S. iMengel et al. 
(I2OO2I ) measure a dynamical mass of 3.3 x 10^ M0, which suggests a steeper than Salpeter 
IMF slope (2.5 rather than 2.35). The discrepancy between ages determined from Br7 
and those determined from photospheric features of supergiants cannot be explained by our 
different IMF choices or different magnitudes or extinctions (these are the same). It may 
be due to different values of EW[Br7] determined from different types of data: a spectrum 
in this work versus narrow-band imaging in theirs, which are subject to different difficulties 
in removing diffuse stellar and nebular backgrounds. The discrepancy may also be due to 
diverging model predictions for our differing spectral diagnostics. Strong Br7 from O stars 
and CO absorption in supergiants are not found simultaneously in the spectrum of a single- 
aged stellar population based on the latest tracks used in StarburstQQ. However, these are 
single-star evolutionary tracks, so the models do not yet account for binary systems (in 
which most stars are thought to reside) or the interactions between binaries and stars in 
cluster environments, e.g. mass transfer, wind-wind interactions, etc. Finally, SSCs may 
not form instantaneously, but rather over the course of several Myr, and if SSCs observed at 
large distances are actually unresolved clusters of clusters then they could have larger age 
spreads-on scales similar to the ELC sizes, 30 Doradus has an age range of about 25 Myr 
due to propagating star formation: R136 c ontains low-mass stars w ith ages of 4 — 5 Myr 
and massive stars with ages of 1 — 2 Myr (IMassey fc Hunteii Il998bl ), while nearby Hodg e 
301 (~ 40 pc away) harbors supergiants with ages of 20 — 25 Myr (e.g. iGrebel &: Chull2000l ). 
We view ELC F and 30 Dor as a warning that ages determined via single-star, single-aged 
population synthesis models may have systematic errors of a few Myr at ages near the end 
of the O star era, and proceed with this caveat in mind. 



4. Mass Loss from ELCs 

With mass estimates for the ELCs in hand, we now assess the ability of the ELCs to 
gravitationally bind their H ll gas by comparing the observed Br7 line widths anii of ELCs 
with their escape velocities Wesc, which depend upon cluster mass and size. We can write Vesc 
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in terms of the ID stellar velocity dispersion of a virialized cluster as Vesc = V^cr*- The 
virial relation can be expressed as al = GM/rjrhp, where rhp is the projected h alf-light radiu s 



and T] is a constant that depends upon the mass distribution of the cluster (e.g. lSpitzerlll987l ). 



The value of t] is often taken to be ab out 10, although for a range of realistic globular cluster 



concentration parameters, 0.5 — 2.5, iMengel et a. 



concentrated clusters have smaller rj). iBoily et al. 



found that t] = 9.7 — 5.6 



more 



(l2005r ) calculated the time evolution of r] 



in SSCs using a gas-dynamical model that incorporates stellar evolution; they found that in 
low-density clusters, r) does not change much over the first 10 Myr, but in intermediate- and 
high-density clusters that resemble M82-F and R136 in 30 Doradus, respectively, rj increases 
by factors of about 1.3 and 1.7, while the clusters grow more compact and r^p decreases. 
An increasing t] implies a decreasing Vesc, which makes it progressively easier for H ll gas to 
escape from a cluster. If the measured anii equals or exceeds Vesc, then the cluster cannot 
gravitationally bind the high- velocity ionized gas in the line wings, which has sufficient energy 
to escape from the cluster potential in an unbound outflow (in the absence of a confining 



extern al pressure). For the two ELCs (S and F) whose cr* were measured by iMengel et al. 



(120021 ). we can directly compare Vesc and anii- Assuming that the measured represents a 
virial velocity dispersion, ELC S has Vesc = 20.1 km s~^ and cthii = 31.8 km s~^ (corrected 
for the instrumental line profile. Table [3l), so most of its gas is not bound to the cluster. If 
the other ELCs in Table [3] share this typical value for an SSC stellar velocity dispersion, then 
nearly all of them have predominantly unbound gas. ELC F has the second-lowest intrinsic 
nebular line width at anii = 11-6 km s~^, and it has Vesc = 28.6 km s~^, which suggests 
that most of its gas is bound to the cluster. If mass segregation is present in these ELCs, 
then the measured cr* would be an underestimate of the true value because it is dominated 
by the light of massive, centrally concentrated stars. In that case, ELCs S and F would be 
better able to bind their nebular gas. 

In order to estimate escape velocities for the other ELCs, we consider the mass of both 
the cluster's stellar component (implied by Starburst99 for a Kroupa IMF extending from 
0.1 to 100 Mq, § [3]) and its ionized gas component (whose mass ranges from 0.01 to 2.6 
times the stellar mas^). We adopt an extinction correction where possible (as in Table H]), 
and we assume that stars and gas fill the same volume out to the smaller of the half-light 
radii measured from the 2.17 and 2.25 /im FWHMs (see Figure |6] and Table [1]). These 
assumptions are conservative in the sense that they place the total amount of mass available 
in the minimum volume, which leads to an upper limit for Vesc- Adopting a value of rj = 10, 
we derive escape velocities for most sources that are well below the observed H ll line widths: 
Vesc/o'Hii ranges from 0.1 to 0.3 for all ELCs except the two clusters with the narrowest Br7 



■^The ionized gas mass is deduced from outflow model fits in Paper II, §2.3 and Table 1. 
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emission (i.e., all bELCs have anii > Vesc)- The exceptions are ELC F, which we have 
already concluded can bind its gas, and ELC Dl, which has the narrowest Br7 line in the 
sample, and Vesc/o'Hii = 2.1. A more conservative estimate of the gas boundedness comes 
from assuming that most of the ELC mass is actually concentrated in a 4 pc radius (median 
value for SSC effective radii in IWhitmore et al.l 119991 ) that is unresolved in our data; this 
leads to higher values of Vesc/crHii = 0.2 — 0.8 for bELCs, but still suggests that much of the 
HII gas in the observed profiles is not bound to the clusters. If we consider both this smaller 
radius an d decrease r) to the m inimum value of 5.6 that corresponds to a high-concentration 
globular (IMengel et al.l |2002| ) . then the ratios Vesc/o'Hii for ELCs Bl, B, and E5 rise to 
0.9 — 1.1, while those of the remaining bELCs are in the range 0.2 — 0.8. Thus even with the 
most compact mass distribution that we imagine they could have, even in the ELCs with 
the largest Vesc/ o'hii ratios, the high- velocity gas in the broad line wings would be unbound. 
This suggests that a significant fraction of the ionized gas in bELCs is unbound and may 
be outflowing, but that narrow-lined ELCs can gravitationally bind their HIIRs and do not 
drive outflows. 

An exception t o this rule may be the compact, embedded broad-line source in NGC 5253 
(ITurner et al.l 120031 ) that resembles the Antennae bELCs in terms of line width, mass, and 
luminosity. Although its Br7 line is very broad (76 km s~^ FWHM) and hence supersonic, it 
appears to be compact enough to bind its H ll gas if the most conservative choices for t] and 
half-light radius (th e smallest radio semimajor axis is 0.5 pc) are adopted for a 10^ Mq cluster 
(ITurner et al.ll2003l ): in this case, Vesc = 131 km s~^, which greatly exceeds the measured line 
width, and leads to the conclusion that the cluster easily binds its UDHII. However, if we 
assume ?7 = 10 and use the .ft^-band size of the nebu la to delineate the Br7-emitting region (1 
pc radius for their H2, lAlonso-Herrero et al.ll2004l ) rather than the radio size, then we find 
Vesc = 69 km s~^ This would suggest that plenty of gas in the Br7 line wings has adequate 
energy to escape the cluster's potential well. Either way, the NGC 5253 cluster is far more 
compact than the Antennae ELCs; this difference may be due to a higher ISM pressure in 
NGC 5253 that confines the UDHII at high density and prevents any significant outflow of 
gas. 

In the foregoing discussion we have assumed that the stars in ELCs are distributed 
like their infrared light. The i^'-band half-light radii of ELCs range from about 10 to 50 
pc, measured either in the continuum narrow-band filter at 2.25 fim, which is dominated 
by starlight (Tabled]), or in the 2.17 fim filter, which also includes Br7 nebular emission 
(Figure E])- ELC radii are up to an order of magnitude larger than the median effective radius 
of optically visible SS Cs of all ages in t he Antennae (4 pc), measured from HST images 
(IWhitmore et al.lll999l . and see Figured]). iMengel et al.l (120021 ) measured an effective radius 
of 3.6 pc for the apparent optical counterpart of ELC F from the /-band HST observations. 
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and they found simi lar values up to 6 pc for a few other optically revealed clusters. However, 
Mengel et al. J2005h reported that the mean /-band effective radius is 16 pc (with a standard 



deviation of 15 pc) for seven isolated young (age < 4 Myr) Antennae clusters. The 2.25 /xm 
half-light radii of ELCs Bl, B, and D, whose ages are below 4 Myr, are 29, 49, and 40 
pc, respectively (Tabled], Figure [6]). Thus young clusters are intrinsically large at visible 
wavelengths even at the WFPC2 resolution of O'.'l. Most ELCs are resolved in our near- 
infrared images (with point source FWHM of 0'.'4 in the combined frames, see Fig. [5]), but 
what we view as a single source could break up into multiple sources when viewed with better 
angular resolution. In order to check this, we examined archival optical and near-IR HST 
images of several ELCs together with ISAAC ones: Figure [2] shows that while some ELCs 
consist of multiple sources, even at the finest spatial resolution they tend to be dominated 
by a single bright one. ELC B, for example, is comprised of a very bright source in K 
band surrounded by fainter emission and at least one neighboring source to the southwest, 
but there are a handfu l of spatially coincident sources in HST filter F555W (Figured] inset, 
Whitmore et al.lll999l ): it may be a cluster of clusters. ELC Bl is a single bright source 
(like most other ELCs) in K band and a very faint one in /. Caution must be exercised 
in attempting to match sources observed at different wavelengths and resolutions, especially 
when extinction can be important. 

The position-velocity diagrams in Figure H] show that the Br7 emission in ELCs Bl, 
B, and D has spatial and velocity structure. ELC C, which appears complex and diffuse in 
optical and near-IR images, also has extremely extended spatial and velocity structure in 
Br7 with a very weak continuum that is centered between its high- and low- velocity 'lobes'; 
it may comprise several cluster components, or it may be driving a bipolar outflow. Figure [6] 
shows that most (but not all) bELCs have slightly more extended emission in the 2.17 /im 
filter, which includes Br7, than in the 2.25 /xm filter, and that there is a weak tendency 
for the ratio of 2.17 /xm to 2.25 /im sizes to increase with age. An HIIR might have 
more extended line than continuum emission (which is dominated by stars) because Lyman 
continuum photons are escaping to ionize the surrounding ISM, and this may produce some 
of the observed extended emission of ELCs. However, another process that would produce 
extended line emission is the escape of ionized gas in an outflow, which is suggested by 
the high velocities of some of the extended Br7 relative to the line centers, as well as the 
double-lobed appearance of ELC C in Figure |U An increase in Br7-to-continuum size ratio 
with bELC age could be due to expansion of an outflow into the surrounding ISM and/or 
contraction of the stellar cluster with time (see § |6]). 
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The Nature of ELC HIIRs 



It is useful to compare the HIIRs that are excited by ELCs with three familiar catagories 
of HIIRs (Galactic compact and ultra-compact HIIRs ((U)CHIIs), Galac tic HIIRs, and extra- 
galac tic GHIIRs), as well as with the recently discovered class of UDHIIs ( iKobulnicky fc Johnson 
19991 ). We will discuss four properties that distinguish different types of HIIRs: luminosity 
(in the form of Q[H^]), line width, size, and electron density He- 

We presented Q[H^] values for ELCs in §[3] and discussed the Br7 line widths of ELCs 
in § HJ we now place bELCs (excluding the narrow-lined ELCs Dl and F) in context with 
other HIIRs by plotting their Q[H'^] vs. H I recombination line widths (with instrumental 
and thermal components, assuming a 10^ K HIIR, removed), together wit h data for other 
objects in Figure [71 The luminosity of 30 Doradus in the LMC is taken from IWalborru (jl99ll ) 
an d its integrated H a line width inferred from the spatially resolved echelle spectroscopy 

(119941). The lurn inosity and Br7 width for the ELC-like source in 
J2003h. For the UDHI Is in H e 2-10 the Q[H+] values are 



of Chu fc Kennicut 



Turner et al. 



NGC 5253 are from 

derived from the radio data of I Johnson fc Kobulnicky I (120031 ) . and the line width adopted 
is the typical value of about 55 km s~^ observed for Br7 (corrected for instrumental and 
thermal contributions, which resembles bELC line widths) in several slits that overlap with 
sever al of the radio so u rces (IHenry et al.l 120071 ). The line widths for a sample of (U)CHIIs 
from iGaray &: Lizand (119991 ) are for H I radio recombination lin es. The sample of top- 
ranked GHIIRs in nearby galaxies is from lArsenault &: RoyI (Il988l ). who compiled Ha line 
widths, luminosities, rms ng values, and diameters from the literature. We assume case B 
recombination to determine Q[H~^] values from line fluxes. Figure [7] shows bELCs at the 
high end of the range in Q[H~^] and FW HM, where they overlap w ith the GHIIRs, which are 
characterized by supersonic linewidths (jSmith &: Weedmanlll970l ) and Q[H~^] ~ 10^° — 10^^ 
Galactic (U)CHIIs have much lower lumin osities because they ar e powered by just one 



or a few massive stars {Q[H^] ~ 

1044_io50 s-i. lCarav fc LizanolEgggI ). but their line widths 



overlap with those of the brighter HIIRs. Not shown are common Galactic HIIRs, which 
exhi bit narrower H recombi nation line FWHMs of 20 — 40 km s"^ and Q[H^] = 4 x 10^^ — 10^^ 
s~^ (IChurchwell et al.lll978l ). A loose correlation between FWHM and Q[H~^] is evident, and 
its slope for bELCs is consisten t with that reported for GHIIR populations of nearby spirals 
(i.e. NGC 4449 and NGC 4321 . lArsenault et al.lll99ol : lRozas et al.lll998l : lFuentes-Masip et aL 



20001). 



We also compare measurements of size (diameter d) and Ue for several types of HIIRs 
in Figure [3 The variations in d and Ue within the (U)CHII, UDHII, and bELC populations 
extend along tracks of constant Q[H~^] for a constant-density Stromgren sphere: Ue oc d~^'^ 
(solid lines in Figure [7]), while GHIIRs cover a broader range in Q[H~^] values. While 
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GHIIRs are similar to bELCs in terms of luminosity or inferred stellar mass, they are larger 
((j^lO O pc) and much les s dense (e.g. ELC Bl has rif. = 10^'^ — 10^ cm~^ from [Fe ill] 



ratios, [Gilbert et al.l 120001 . and see below). GHIIRs have rms densities of a few to a few 
100 cm~^ (derived from emission measure profiles), but for bidden-line diagri ostics, e.g. [S ll] 
and [O ll], yield densities that are ten times greater (e.g., iKennicuttI 1 1984 ) . This range of 
densities is commonly interpreted as evi dence for dumpi ness of the HIIRs, and it suggests 
filling factors of dense gas of 0.01 — 0.1 (IKennicuttI Il984j ) . Such a variable-density medium 
i n GHIIRs is certainly p lausible: inhomogeneities are common in resolved HIIRs like 30 Dor 
(IChu &: Kennicuttlll994j . e.g.), and density variations have been me asured in both Galactic 
and extragalactic HIIRs (ICopetti et al.ll2000l : ICastaneda et al.lll992l ). 



(U)CHIIs are included in Figure [7] because they represent extremes in all of the plotted 
parameters. Although their low Q\H^\ values, small implied stellar masses, and sub-parsec 
sizes clearly differentiate Galactic (U)CHIIs from GHIIRs and bELCs, their extremely high 
densities are similar to that of ELC Bl: compact (diameter d < 0.5 pc) and ultra-compac t 



{d < 0.05 pc) HIIRs have Ue ranging over ~ 10^ — 10^ cm ^ (jWood fc Churchwelll Il989l ). 



Common Galactic HIIRs (not plotted) have comparable luminosities to those of (U)CHIIs, 
but they are larger (0.5 — 15 pc) and have slightly lower densities (10^ — 10^ cm~^) than 
their compa ct cousins. Moreover, UC HIIs have dynamical ages indicating extreme youth 
(5000 years. IWood fc Churchwelll Il989l ). although they could live much longer (~ 10^ years) 
(iGaray fc Lizandll999l ). so the UCHIIs may be the precursors of common HIIRs. The bELCs 
appear to be 3 — 7 Myr old (§ [3]) and therefore much more evolved than UCHIIs, yet they are 
clearly more compact and perhaps at an earlier stage of evolution than the diffuse GHIIRs. 

We estimate a volume-averaged ne for all ELCs by assuming ionization and recombi- 
nation balance for a 10^ K HIIR with Q[i/^] (extinction-corrected, where possible) from 
Table m along with their Br7 sizes (i.e. the FWHM measured in our 2.17 /xm narrow-band 
image, see Figure E]). The resulting values range from 30 to 600 cm~'^ with an average of 
180 cm~^, and they represent lower limits on the electron deii s ities o f ELCs. For the four 
bELCs that have thermal radio counterparts, iNeff fc UlvestadI (120001 ) compute similar rms 
Ue values of 30 — 360 cm~^ by assuming measured diameters of up to 70 — 120 pc (which in 
some cases gives an upper limit o n size and lower l imit on Ue because the sources are not 
resolved by their ~ 80 pc beam). iNeff &: UlvestadI (120001 ) also consider the possibility that 
their thermal radio sources have 8 — 30 pc diameters that are nearer to those observed for 
SSCs (~ 8 pc. lWhitmore et al.lll999l ). and thereby obtain higher estimates of ne = 1 — 5 x 10'^ 



cm~^. If the i^-band emission from ELCs is concentrated in compact, unresolved cores that 
are the size of the optical SSCs, then the density estimates would increase correspondingly, 
to 700 — 8500 cm~^ with a mean of 2250 cm~^; the maximum value is that of ELC Bl, and 



it is consistent with the range inferred from [Fe ill] ratios by lGilbert et al.l (120001 ). Thus the 
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radio-inferred densities for a subset of ELCs are comparable with the IR-inferred ones. For 
consistency, the ng values plotted in Figure [7] are the lowest (rms ne) estimates for GHIIRs, 
ELCs, and UDHIIs, which we discuss next. 

We estimate densities for the UDHIIs o f He 2 — 10 in the same way a s for the ELCs, 
using t he Q\H~^] values fro m 0.7 cm data of I Johnson fc Kobulnicky I (120031 ) and the 10 /xm 
sizes of IVacca et al.l (120021 ). where we combine their objects 1 and 2 since they are blended 
in the observations of the latter work. Thus we assume that the mid-IR and near-IR sizes 
are comparable, and that the stars and warm dust share the same volume. We find densities 
for the two compact UDHIIs (4 and 5, with radii of 4 and 5 pc) of 2200 and 1280 cm~^. For 
the two larger UDHII regions (1+2 and 3, with radii of 16 and 17 pc), we find Ue ~ 390 and 
180 cm~^. The lowest estimates of average ELC densities are thus similar to those of the 
more extended UDHIIs in He 2 — 10, and the SSC-size estimates are similar to those of the 
compact UDHIIs. 



However, IVacca et al.l (120021 ) do not assume that the stars and warm dust in UDHIIs 
share t he same volu r ne; th ey estimate rie for the He 2 — 10 UDHIIs by fitting the radio-NIR 
SEDs. IVacca et al.l (120021 ) model a UDHII as a central HIIR that envelopes a young SSC 
and extends out to the inner radius of a surrounding annular dust cocoon that extends to an 
outer radius. Their best-fit inner radii are similar to the siz es of optical SSCs, 2 — 5 pc, an d 
they are consistent with fits to more recent radio data of I Johnson fc Kobulnicky I (120031 ). 
This leads to average values for of 1400 — 5720 cm~^, wh ich a gree well with the radio 



estimates of the average from iKobulnicky fc Johnson] (119991 ) and jjohnson fc Kobulnicky 



Turner et al 



J2000h 



(120031 ). For the compact radio counterpart of the NGC 5253 HIIR, 
report an even larger estimate of the rms ne of 4 x 10'^ cm~^, which is in the range of UDHII 
Tie. Adopting the higher estimates of for UDHIIs or bELC^ would move them toward to 
the (U)CHIIs in the Ue — d panel of Figure [71 while adopting the higher estimates for GHIIRs 
would move them up; the figure would still display a strong anticorrelation between Ug and 
diameter in which bELCs and UDHIIs occupy a central region that is distinct from that of 
the other massive systems plotted, the GHIIRs. 

Although the populations of HIIRs that are discussed above exist in a variety of envi- 
ronments and were all observed in different ways and subject to different selection effects 
from ELCs, their comparison shows how ELCs fit into the broader context of star formation. 
Their nearest relatives are the UDHIIs, whose Lyman continuum luminosities and electron 



''We discuss further evidence for high-electron-density material in bELCs in Paper II, where we derive 
radial density distributions for the bELC HIIRs by fitting the Br7 profiles with a kinematic model for a 
cluster outflow; this yields even higher values for the emissivity- weighted average Ue- 10^ — 10^ cm~^. 
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densities are similar, but UDHIIs are so heavily obscured that they are rarely visible in the 
optical range and are thought to be younger than ELCs. UDHIIs are classified based on their 
radio properties, and ELCs based on their recombin ation lines, b u t we h ave already confiated 
the two groups by assuming the Br7 line fluxes o f iHenrv et al.l (120071) ern erge from the ra- 
dio/MIR sources of iKobulnicky fc JohnsonI (119991 ) and IVacca et al.l (120021 ) . That is because 
all evidence suggests that ELCs and UDHIIs represent the same objects, perhaps at slightly 
different evolutionary stages, so that massive clusters form within gas- and dust-shrouded 
UDHIIs, which become ELCs as their massive hot stars drive winds (that are visible in re- 
combination line emission) and expel gas to reveal young SSCs with low extinctions at visible 
wavelengths. If t hese massive cluster s form hierarchically from clusters of subclusters (e.g. 
Testi et al.ll2000l : iBonnell et al.ll2003l ). then the ELCs may be comprised of several evolving 
UDHIIs that our data cannot resolve. We conclude that bELCs and UDHIIs have some 
properties in common with other types of HIIRs, but taken together, their broad lines, high 
densities, intermediate sizes, and large luminosities constitute a new class of massive HIIRs 
that is distinct from GHIIRs. 

While our small sample spans two orders of magnitude in Q[H'^] and inferred clus- 
ter mass, it may represent only the bright end of the full young cluster population in the 
Antennae. If the mass function of young clusters in the Antennae starburst reaches down 
to single-star scales, would those objects resemble Galactic UCHIIs? The analogy between 
UCHIIs and UDHIIs may extend beyond extreme densities, compactness, and extinctions: 
the massive stars within UCHIIs are thought to spend 10 — 20 % of their main-sequence 
lifetimes in an embedded phase ( iWood &: Churchwelll Il989l ). and those in UDHIIs may 
spend even longer fractions o f their lives hidden from optical view, e.g. 15 % in He 2 — 10 
(iKobulnickv fc JohnsonI 1 19991 ) to 40 % in NGC 5253 jMartm- Hernandez et allbood ). until 
their cluster winds remove much of the obscuring gas. Such an extended embedded phase 
for massive stars in starbursts is proposed to explain discrepancies between measured mid- 
IR line ratios a nd st ellar population synthesis models with normal upper mass limits by 
Rigby fc Riekd (120041 ). who point out that the higher ISM pressures in starbursts should 
lengthen their embedded compact HIIR phase relative to that in the Galaxy. The impli- 
cations of a long embedded phase in whi ch IR nebular lines a re suppressed are serious for 
starburst population synthesis modeling (iRigby fc Riekd l2004j ). 



6. Evolution of ELCs 



In an evolutionary sequence determined by cluster age and embeddedness, massive stel- 
lar clusters form in heavily obscured UDHIIs, evacuate their surroundings as wind-blowing 
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ELCs, and, if they survive the mass-loss phase, emerge as SSCs, which may ultimately 
become GCs if they survive further destructive cluster dynamical processes and tidal in- 
teractions with their host galaxy. This scenario requires some explanation of cluster size 
evolution, in particular how the large K-hand radii of ELCs (10 — 50 pc. Table [T] and Fig- 
ure [6]) can be r econciled with the typi cal half-light radii of optically revealed SSCs (~ 4 pc in 
the Antennae, IWhitmore et al.lll999l ). We assume that ELCs are young SSCs embedded in 
HIIRs because their inferred stellar masses are similar to those of SSCs, which can be factors 
of a few to ten times more compact than ELCs. FigureOshows a wide scatter in the 2.25 fim 
stellar continuum sizes of 10 ELCs at ages of 6 — 7 Myr: their mean (± standard deviation) 
half-light radius is 17 ± 8 pc. Younger ELCs (5 have ages of 3 — 6 Myr) have a larger mean 
half-light radius of 35 ± 11 pc at 2.25 fim. Additional evidence for a decrease in Antenna e 
cluster size with age, as well as a large scatter in size, is presented by iMengel et al.l (120051 ) . 
who report average (=b standard deviation) WFPC2 /-band effective radii for 18 isolated 
clusters with ages below 4 Myr and in the range 8 — 11 Myr of 16 ± 15 pc and 6.5 ± 5.3 pc, 
respectively. 

This tendency of older clusters to be more compact could be explained by cluster con- 
traction with age or by the preferential survival of more concentrated clusters at any given 
age. Clusters that form hierarchically would contract as subclusters merge to form a single 
cluster. Simulations by lBonnell et al.l (120031 ) show that star clusters may form hierarchically, 
such that fragmentation in a molecular cloud produces sheets and filaments which form stars 
before the rest of the cloud material, and these stars form high-density subclusters that 
merge after several free-fall times to create a centrally condensed cluster. Observations of 
the S erpens molecular cloud core provided the first evidence for hierarchical cluster forma- 
tion: iTesti et al.l (120001 ) find that star formation is not homogeneous in the core, but rather 



is concentrated in dense subclusters that are well separated spatially and kinematically. If 
ELCs are young enough to be undergoing the final merger of their initial subclusters, this 
might explain their size trend in Figure O 

However, even if hierarchical formation is still underway in ELCs and causing contrac- 
tion, it must compete with the countervailing influence of mass loss due to stellar winds that 
lead to cluster outflows. Cluster expansion is caused by the loss of gas via stellar evolution 
and stars via relaxation, tidal stripping, and dynamical friction. Mass loss and consequent 
expansion can unbind a cluste r and cause it to dis solve below detection limits into the back- 
ground stellar population (see lFall fc Zhangll200ll . for an analysis of these effects across the 
cluster mass spectrum); hence compact clusters live longer than less compact ones, and they 
increasingly dominate the cluster population as they age. Two recent studies of the age 
distribution of stellar clusters in the Antennae argue that the sharp decrease in the number 
of clusters as a function of age is strong evidence for rapid dissolution of the majority of 
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young clusters on timescales of order 10 Myr (IMengel et al.ll2005l : iFall et al.l 120051 ) . Thus 
we favor expansion and rapid dissolution over contraction due to hierarchical formation in 
explaining the smaller sizes of older ELCs. 

Large ELCs may therefore be puffed-up young SSCs that have expanded due to prodi- 
gious mass loss, which limits their chances for survival if they are not compact enough, 
if their IMF is top-heavy or truncated at l ow rn asses, or if their star-formation efficiency 
(SEE) is not high enough. iBoily fc Kroupal (120031 ) simulate cluster evolution in the presence 
of rapid (i.e. instantaneous) mass loss; after equilibrium is re-established, they find expan- 
sion factors for the half-mass radii of up to a factor of ten or more for bound clusters having 
Plummer and King distribution functions with varied concentrations and assumed SEEs. 
This is a plausible explanation for the sizes of ELCs if they were born as more compact 
clusters within UDHIIs; however, the mass loss in ELCs proceeds over the course of at least 
several Myr, so t hey w ould not suffer such extreme expansion as the clusters simulated by 
Boilv fc Kroupal J2003h . 



The expansion of young clusters due to dynamical evolution is observed in nearby re- 
solved clusters: cluster core r adii in the LMC and SMC increase with age, at varying rates, 
up to ages of at least 1 Gyr (lElson et al.lll989l : iMackey fc Gilmorell2003a| ]bl). These young 
clusters (10® — 10^ yrs old) have power-law luminosity profiles that often show bumps or 
dips that suggest the presence of substructures that have not been smoothed out by re- 
laxation ( Elson et al. 1987 : Mackey fc Gilmore 2003bl lal). These profiles do not display the 
tidal truncation of a King profile ( King! 1966 ). which best describes the luminosity profiles 
of much more evolved GCs. There is evidence that the more massive SSCs in nearby galax- 
ies also have power-law light profiles, some extending beyond 100 pc in radius, some with 
clumpy substructure (see references in ISchweizerl 12004 ) : evidence for the onset of tidal trun- 
cation in older SSCs (beyond a few 100 Myr) and interm ediate-age GCs has also been found 
(IWhitmore et al.lll999l : ISchweizerll2004l ). iLarsenl (120041 ) finds that extended envelopes rela- 
tive to King profiles are common among the youngest clusters in a sample of spiral galaxies 
observed with HST WEPC2, but his sample does not show the trend of increasing cluster 
radii with age that is seen in the LMC and SMC. 

The interpretation of large ELCs as puffed-up SSCs that are caught in the act of dis- 
solving assumes that their K-hand light at large radii is dominated by stellar continuum 
emission. Alternative explanations for their sizes might, however, involve dust scattering 
and emission from a shell at larger radii surrounding an SSC-sized cluster, or differential 
extinction that preferentially obscures the cluster core relative to its outer radii, making the 
observed EWHM larger than the intrinsic one. A final possibility is that SSCs may form hi- 
erarchically from clusters of smaller subclusters that merge rapidly to make a single massive 
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cluster (within a couple of free-fall times for a low-mass cluster, iBonnell et al.ll2003l ): in this 
picture, prior to the final merger into a single cluster, the average stellar density is smaller 
than in individual subclusters, and the overall structure is probably not spherically symmet- 
ric. Observations at higher angular resolution at infrared wavelengths, either space-based or 
adaptive optics-assisted, will be able to test these scenarios. 



7. Star Formation in Massive Clusters 

Comparing the star-formation rate (SFR) in ELCs with that of the entire system is 
difficult because the SFRs of galaxies are determined from observed quantities such as Ha 
and far-infrared (FIR) luminos ities averaged oy er large spatial scales, and the y consider 



contin uous star formation (e.g. iKennicuttI 119981 ). Assuming the calibrations of iKennicutt 



( 119981 ). ad opting an IRAS lum inosity of 10^^ for the Antennae leads to an SFRfir ~ 17 



Mq yr-^ 



Zhang et al.l (120011 ) use Kennicutt's Ha calibration to estimate SFRg 



a 



20 



Iq yr If we apply this Ha calibration to the sample of ELCs by assuming case B 
recombination, we derive a total SFRj^^ of 3.2 M© yr~^. Thus the fraction of the recent 
SFR of the Antennae that is provided by ELCs is at least 15%, and half of this contribution 
is provided by ELC Bl alone. This is a lower limit because we have neglected extinction 
in two-thirds of the ELCs, our sample probably does not include all ELCs, and there could 
be a significant escape fraction for ionizing photons from the clusters. From a ph otometric 
study of the current star cluster population in the Antennae, iMengel et al.l (120051 ) estimate 
a total SFR in clusters over the past 25 Myr of at least 20 M© yr~^ (in agreement with the 
Ha estimate for all star formation), but they argue that the age distribution of SSCs implies 
rapid dissolution of most clusters, and that accounting for this could raise the SFR to at 
least 50 M© yr'^ 

An estimate of the fraction of stars that are formed in clusters rather than in the field 
may be derived from imaging observ ations that conipare amounts of discrete and diffuse 



emission from star-formation tracers: iFischer et al.l (119961 ) reported that about 70% of the 



Br7 emission in a 30" aperture including the nucleus of NGC 4039 and the southern overlap 
region is diffuse, while the remain ing 30% is associated with Ha and radio knots. Similarly, 



Hummel fc Van Per HulstI (119861 ) report that the total fraction of radio continuum emission 
from discrete knots (at 6" resolution) is about 35%, with errors of 20 — 30%. Our lower limit 
to the fraction of the star-formation rate in ELCs is about half of what one might expect 
from these imaging studies. These results are consistent with the range deduced from UV 



with stellar clusters (e.g. 



Meurer et al. 


1995; 


Maoz et al. 


1996) 
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of the Wolf -Rayet starburst He 2-10 rev eals that 15% of the O stars in the galaxy reside 
in UDHIIs ( Kobulnicky &: JohnsonI [19991 ). which are probably the progenitors of ELCs and 
SSCs. Since our measurement is a lower limit on the fraction of star formation in Antennae 
clusters, it appears that ELCs constitute the majority of the clustered star-formation rate in 
the Antennae, and at least 15% of the Ha-derived total SFR. 

The actual proportion of clustered star formation in the Antennae may be nearer to 



100% if many of the y oung SSCs are quickly dissolved or disrupted (e.g. iMengel et al. 



20051 : iFall et al.l l2005l ). IWhitmore et al.l (119991 ) find a power-law form for the luminosity 



function for Antennae SSCs that c 



(see e.g. 


Meurer 


1995: 


Larsen 


2002 


) as well as 


Kennicutt et al. 


19891: 


McKee & Williams 


1997) 



and o bservatio ns have begun to re veal t he process fe.g. Ide Grijs et al. 



[e.g. 



ve into a log-norma l 



into the field population in order for a power-law SSC population to evo 
globular cluster population; simulations indicate that thi s is plausible (Fall fc Zhanej 120011). 



2003 



Lamers et al. 



20051 ). Indeed, iMengel et al.l (120051 ) and iFall et al.l (120051 ) argue that rapid dissolution into 
the field population is the normal fate for SSCs in the Antennae; those that survive mass loss 
in the ELC phase are subject to further destructive dynamical and evolutionary processes 
that may prevent them from becoming GCs. 



8. Conclusions 

High-resolution spectroscopy of the youngest SSCs in the Antennae reveals broad, spa- 
tially extended, non-Gaussian, supersonic Br7 emission (60 — 110 km s~^ FWHM) in the 
ELCs. Typical GHIlRs also have supersonic linewidths and luminosities at the low end of 
the ELC ranges, but they are much larger and less dense than ELCs, and are generally ex- 
cited by lower-mass clusters and OB associations. Thus ELCs appear to represent a different 
class of dense extragalactic HllR that is powered by young SSCs and may be associated with 
high-pressure environments. Their closest relatives are the UDHIIs, which are younger (age 
^1 Myr), more heavily embedded radio and mid-IR nebulae with high electron densities. 
UDHII Lyman continuum luminosities and infrared sizes are similar to those of ELCs, so 
UDHIIs require similar stellar masses for their embedded ionizing clusters as ELCs, 10^ — 10^ 
Mq. Thus UDHIIs are the probable precursors of ELC. 

The youngest ELCs (ages 3 — 6 Myr) are concentrated in the southern overlap region 
of the Antennae, while those in the northern overlap region are slightly older (6 — 7 Myr). 
The rms radial velocity dispersion of the youngest ELCs is consistent with that measured 
for molecular gas via CO in the same region, so the young clusters and dense gas share 
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similar kinematics. ELCs in our incomplete sample comprise about 15% of the current SFR 
in the Antennae (half of that is concentrated in the brightest cluster) and they make up a 
significant fraction of the expected clustered star formation in the system. 

Because bELC Br7 hnewidths typically exceed the most conservative estimates of the 
cluster escape velocities, much of the H ll gas is not bound to the clusters-it is flowing out 
in a chistcr wind. The greater extent of Br7 emission relative to continuum emission is 
consistent with escape of ionized gas in winds. Although this mass loss takes place over the 
course of several Myr rather than instantaneously, it may be dramatic enough to prevent 
some of the ELCs from surviving as bound SSCs; their large A'-band sizes and the decrease 
in mean ELC size with age may be the signature of rapid cluster expansion and dissolution, 
which leads to the preferential survival of the most compact clusters in each age group. In 
Paper II we present a kinematic model for cluster outflows to infer their mass-loss rates and 
energetics, and discuss ELC feedback in the context of similar systems such as wind-blown 
superbubbles and galactic superwinds. 
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Fig. 1. — Image of the Antennae using ISAAC 2.17 fim narrow-band filter, which includes 
Br7. Units are arcseconds of R.A. and Dec; north is up and east is to the left. Diamonds 
mark positions of NIRSPEC targets (except for the circled sources 10 and 16 where they 
are omitted to avoid crowding). Circles mark positions of the brightest young clusters in 
HST F555W (Whitmore et al. 1999). SSC B is shown in an inset, magnified by a factor of 
five, in order to reveal IR structure and separate HST sources. At the origin is Star 4 of 
Whitmore et al.l (|l999l ). 
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Fig. 2. — Comparison of ELCs observed from the ground and by HST in several wavebands. 

Bottom row shows the ISAAC 2.17 image (0'.'4 seeing); subsequent rows show NIC3 
F222M image (O'.'2/pixel), NIC2 F187N and F160W images ((/.'OTS/pixel), and ACS F814W 
image ((X'OS/pixel). Horizontal bar is 1" long. 
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Fig. 3. — Br7 spectra of young Antennae SSCs and nuclear regions, featuring broad, non- 
Gaussian wings. The spectral resolution is about 12 km s~^ FWHM (4 pixels, shown as a 
bar in upper right corner of all panels), so the lines are well resolved and their widths are 
not simply due to thermal broadening. He I emission at 2.1647 /xm (A^; = —191 km s^^from 
Br7) is also visible in the brightest sources (top left panel) . Note the different velocity scale 
in top left panel, and varied flux scales throughout. 
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Fig. 4. — Br7 position- velocity diagrams for ELCs Bl, B, C, and D, overlaid with arbitrarily 
scaled integrated spatial profiles and extracted spectra. The Br7 emission is extended with 
respect to the continuum, and its spatially resolved velocity gradients are suggestive of 
nonspherical flows. 
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Fig. 5. — Normalized, azimuthally averaged radial profiles of ISAAC 2.17 and 2.25 //m 

image counts for two ELCs and Star 4 demonstrate that the ELCs are resolved relative to 
the observed PSF. Error bars indicate la Poisson errors. 
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Fig. 6. — Above: Deconvolved sizes of ELCs measured from 2.25 //m ISAAC narrow-band 
image, plotted as a function of ELC age. At ages of 6 — 7 Myr, the population of ELCs 
has a large size spread (down to half-light radii of 9 pc for ELC F), but younger ELCs are 
larger, with radii of 30 — 50 pc. Below: The ratio of 2.17 /xm (stellar continuum plus Br7 
emission) to 2.25 fim (stellar continuum) half-light radii for bELCs increases slightly with 
age. Most bELCs (which exclude ELCs F and Dl) have larger sizes in the Br7 filter than in 
the continuum one. Sizes are geometric mean FWHMs determined from 2D Lorentzian fits 
with Star 4 geometric mean radius subtracted in quadrature, and error bars are computed 
from fit errors. 
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Fig. 7. — Above: bELC Lyman continuum rate versus line width, with instrumental re- 
sponse and thermal component removed, compared with other samples of HIIRs: nearby 
mass ive cluster 30 Dor (luminosity and line width fr om 



Walborn 



994 respectively), UDHIIs/bELCs in NGC 5253 flTurner et al 



( Kobulnicky fc Johnson! Il999 : 



Johnson &: Kobulnickv 



2003; 



1991 



Chu &: Kennicutt 



Vacca et al. 



2003 ) and He 2 - 10 



20021 : iHenrv et al. 



2007), Galactic (U)CHIIs |Garav &: Lizandll999l ). and first-ranked GHIIRs in several galax- 



ies ( lArsenault fc Roylll988l ). Below: Mean electron density versus diameter for bELCs com- 
pared with the other HIIR populations. Solid lines indicate the Stromgren sphere relation 
for an HIIR with constant density: rie oc d~^'^, for Q[H^] = 10'^^, 10^^, and 10^^ s^^. Sizes 
are measured from radio data for all but the bELCs and UDHIIs, whose sizes are from near- 
and mid-IR images, respectively; adopting the pc-scale radio sizes for UDHIIs or the typical 
optical SSC sizes for bELCs moves them closer to the (U)CHIIs. All densities plotted are 
average values; the GHIIRs would move up by one or two orders of magnitude if we adopted 
the higher densities that are derived from optical forbidden lines and are weighted toward 
dense, presumably low-filling-factor gas. 
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Table 1. Photometry & Sizes of Selected Antennae Sources 



ssc 


A R.A.'^ 


A Dec^ 




FWHM'^ 


Name 


(") 


(") 


(mag) 


(") 


Bl 


-17.7 


-21.0 


14.7 


0.62 


B 


-23.1 


-18.2 


"1 /I 

14.5 


1.05 


D 


-11.7 


-4.5 


15.4 


0.87 


C 


-20.2 


-6.3 


15.2 


r\ AC 

0.46 


Dl 


-13.9 


-2.7 


15.8 


0.23 


D2 


-10.3 


-1.1 


16.9 


0.75 


El 


-11.7 


15.5 


16.1 


0.28 


E2 


-13.0 


14.7 


15.6 


0.54 


E3 


-19.5 


11.4 


17.1 


0.23 


F 


-11.9 


26.6 


15.9 


0.20 


E5 


-8.9 


21.8 


15.3 


0.50 


F2 


-15.9 


32.5 


15.8 


0.36 


Fl 


-13.9 


28.8 


15.7 


0.33 


E4 


-10.7 


20.4 


15.7 


0.77 


Al 


-39.8 


-23.0 


17.5 


0.32^ 


S 


-82.3 


32.1 


14.0 


0.30 



^Positions in arcseconds relative to Star 4 of Whitmore et al. 1999, whose 
position is report ed as R.A. = 12'' 1™ 56.04 ° and DEC = -18°52'43"66 
(J2000), although Iwhitmore fc ZhanJ i2002l) find an offs et of l"2 to the 
southwest of the HST positions from radio observations iNeff &: UlvestadI 
120001 '). All offsets are measured from our NIRSPEC N7 mosaic except that 
of S, which is derived from [Whitmore et all il999| ) . 

magnitude fro m NIRSPEC SCAM photometry, except for SSC S from 
iMengel et all i2002|) . Systematic errors arc 0.2 mag. 

'^Deconvolved FWHM (from geometric mean radius of 2D Lorcntzian fit) 
from distortion-corrected 2.25 /J-m narrow-band ISAAC image; at the An- 
tennae 1 — 93 pc. Lorentzian FWHM of Star 4 is 0.27 , which was used 
to derive seeing- deconvolved FWHMs. 

'^Al is too faint and confused to fit as above, so its FWHM is estimated 
from a spline fit to the azimuthally averaged profile, subtracting in quadra- 
ture Star 4's size measured the same way (0.38 ). 
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Table 2. Log of NIRSPEC Spectra of Antennae Sources 



Source Date texp 

Names (UT) (s) 



B,BP 


2002- 


-02- 


22 


2 


X 


600 


C,D 


2002- 


-02- 


22 


2 


X 


600 


D1,D2 


2002- 


-02- 


22 


1 


X 


600 


E1,E2,E3 


2002- 


-02- 


22 


3 


X 


600 




2002- 


-02- 


23 


2 


X 


600 


W99 10,16= 


2002- 


-02- 


23 


2 


X 


600 


F1,F2,E5 


2002- 


-02- 


23 


3 


X 


600 


E4 


2002- 


-02- 


23 


1 


X 


600 


N4039 nucleus, Al 


2002- 


-02- 


23 


1 


X 


300 


N4038 nucleus 


2002- 


-02- 


23 


1 


X 


300 


S^ 


2000- 


-12- 


11 


7 


X 


900 


W99 r 


2001- 


-02- 


04 


10 


X 


900 



^SSC Bl is SSC A of iGilbert et all fl2000h . 



^ Young cluster 15 of IWhitmore et al 



(Il999f ). 



= Young clusters 1 0, 16, and 1 of 



Whitmore et all fll999f ). 



'^Young cluster 2 of IWhitmore et al.l (Il999l ). 
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Table 3. Br7 Line Measurements 



Name 


F[Br7]xlOi6 ^ 
(erg s"^ cm^^) 


EW[Br7] 
(A) 


b c 

(km s~^) 


Bl 


155.7 ± 2.4 


254.8 ± 8.0 


14/6 


/O 


B 


166.1 ± 1.6 


260.9 ± 5.3 


1447 


63 


D 


47.5 ± 1.7 


189.7 ± 10.6 


1390 


57 


C 


40.5 ± 2.1 


63.8 ± 3.6 


1383 


109 


D1 


3 1+06 


OA 4- 10 A 


1485 


14 




30.9 ± 7.2 


123.2 ± 52.4 


1487 


55 


El 


5.1 ± 0.2 


2.6 ± 0.1 


1599 


40 


E2 


5.2 ± 0.2 


2.5 ± 0.1 


1580 


39 


E3 


5.2 ± 0.3 


1.7 ± 0.1 


1544 


53 


F 


6.7 ± 0.5 


8.4 ± 0.7 


1583 


30 


E5 


6.2 ± 0.5 


20.8 ± 3.6 


1558 


44 


F2 


13.4 ± 1.1 


2.8 ± 0.3 


1601 


80 


Fl 


2.3 ± 0.5 


15.9 ± 3.5 


1566 


42 


E4 


34.5 ± 1.8 


6.0 ± 0.4 


1558 


54 




64.2 ± 1.8 


1.6 ± 0.0 


1575 


45 


N4038 


33.4 ± 3.5 


1.6 ± 0.2 


1597 


77 


S 


6.6 ± 2.2 


4.8 ± 1.6 


1566 


76 


flux 

^Average 
km s~" "'" . 


is not corrected for extinction (see Table [4l . 
error in baryccntric Br7 radial velocity is 1 km s 


^ , with a 


maximum of 2.8 



■^Observed FWHM, not corrected for instrumental line-spread function (FWHM 12 
km s "'"). This correction is significant for the marginally resolved line of SSC Dl. The 
corrected FWHM of SSC F is 27 km s~^. Average error in Br^ FWHM is 2.6 km s^"^, with 
a maximum of 7.5 km s~^. 

'^For D2, total flux and EW arc listed, but Vj-j^^^j and vp^HM refer only to the broad 
component of the line. 

^Refers to an off-nuclear HIIR near the NGC 4039 nucleus. 
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Table 4. Derived Ages & Masses 



IN dlilt: 


T no- DfN+l^ 


Age 


iviabo 






IS ) 




[lU iVlQj 


l^magj 


Bl 


52.71 


3.45 


4.23 


1.2^ 


B 


52.74 


3.44 


5.00 


0.2^ 


D 


52.19 


3.85 


1.93 




C 


52.12 


5.72 


4.13 




Dl 


51.01 


6.13 


1.62 




D2 


52.00 


5.36 


0.80 


0.8"^ 


El 


51.22 


6.75 


0.26 




E2 


51.23 


6.76 


0.41 




E3 


51.23 


6.94 


0.07 




F 


51.34 


6.47 


0.74 


0.1^ 


E5 


51.31 


6.19 


2.60 




F2 


51.64 


6.73 


0.35 




Fl 


50.87 


6.29 


1.47 


0.3^ 


E4 


52.05 


6.55 


0.65 




Al 


52.32 


7.00 


0.50 




N4038 


52.04 


7.00 


3.45 




S 


51.33 


6.61 


3.22 


0.0^ 



^Q[H~^] is derived from Br^ fluxes and not corrected for extinction, which scales 
it up by factor lO^'^'^K . 

'^Ages and masses (uncorrected for extinction) were derived by comparing mea- 
sured K magnitudes and Br'y EWs with Starburst99 models for a cluster with solar 
metallicity and Kroupa IMF (0.1 — 100 Mq). A nonzero Aj^ does not affect ages, 
but does scale masses by factor lO'^ '^^K . 

"^Errors in ages are derived from Starburst99 flts and are typically below 0.1 
Myr. For ELCs D2 and Fl the large errors on EW lead to larger age errors: D2 
has age 5.4^^'^ Myr and Fl has age 6.4^Q Qg Myr. 

'^Relative errors in masses are typically 20% and dominated by photometric 
errors for most sources. Exceptions are ELCs Dl, D2, and Fl, whose relative 
errors reach 37, 30, and 80%, respectively, due to their large age uncertainties. 

° Estimated from Br-/ fluxes and radio measurements reported by 
iNefT &: UlvestadI l |2000l ) (see § [sj . 

^Reported bv lMengel et al.| l |2002| '). 



